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FLUTTER OF A UNIFORM WING WITH AN ARBITRARILY PLACED MASS 
ACCORDING TO A DIFFERENTIAL-EQUATION ANALYSIS 
AND A COMPARISON WITH EXPERIMENT 

By Harry L. Runyan and Charles E. Watkins 


SUMMARY 

A method is presented for the calculation of the flutter speed 
of a uniform wing carrying an arbitrarily placed concentrated 
mass. The method, an extension of recently published work 
by Goland and Luke , involves the solution of the differential 
equations of motion of the wing at flutter speed and therefore 
does not require the assumption of specific normal modes of 
vibration. The order of the flutter determinant to be solved 
by this method depends upon the order of the system of differen- 
tial equations and not upon the number of modes of vibration 
involved. 

The differential equations are solved by operational methods, 
and a brief discussion of operational methods as applied to 
boundary-value problems is included in one of two appendixes. 
A comparison is made with experiment for a wing with a large 
eccentrically mounted weight and good agreement is obtained. 
Sample calculations are presented to illustrate the method; and 
curves of amplitudes of displacement, torque, and shear for a 
particular case are compared with corresponding curves com- 
puted from the first uncoupled normal modes. 

For convenience, the method employs two-dimensional air 
forces and could be extended to apply to uniform wings with 
any number of arbitrarily placed concentrated weights, one of 
which might be considered as a fuselage. The location of such 
masses as engines, fuel tanks, and landing-gear installations 
might be used to advantage in increasing the flutter speed of a 
given wing. 

INTRODUCTION 

The common procedures in flutter analysis of an airplane 
wing involve many simplifying assumptions. In particular, 
the degrees of freedom of the wing are usually determined by 
choosing the first few normal modes of the structure, and the 
wing motion at flutter is then described in terms of these 
chosen modes. This approach of employing prescribed 
modes is often adapted to the Rayleigh type analysis of 
vibration and may be referred to as “Rayleigh type analysis.” 
In specific calculations with this method the amount of work 
required is proportional to the number of normal modes 
involved. In particular, the order of the flutter determinant 
that must be solved depends directly upon the number of 
modes involved. For simple wings, without concentrated 
masses, the Rayleigh type analysis usually yields satisfactory 
results with not more than two or three normal modes. 
However, if the wing carries concentrated masses, such as 


engine, fuel tank, or landing-gear installations, so many 
normal modes may be required to obtain satisfactory results 
that the Rayleigh method may not be the most feasible 
method. 

In cases where many degrees of freedom are involved the 
most logical procedure would be to treat the system of 
differential equations of motion of the wing rather than to 
choose specific modes. For arbitrary wing plan forms this 
method would be in general very difficult and tedious to 
carry through, although it has the advantage that the order 
of the flutter determinant that must be solved depends only 
upon the order of the system of differential equations and not 
upon the number of modes of vibration involved. 

As early as 1929 Kiissner (reference 1) used the differential 
equation approach to formulate the problem in the form of 
an integro-differential equation for a wing of general plan 
form. Kiissner set up some particular examples and sug- 
gested a method of solution by a process of iteration. This 
method was not followed up until during the war when some 
related work was undertaken in Germany but not finished. 
Wielandt (reference 2) has recently made contributions to the 
treatment of nonself-adjoint differential equations by itera- 
tive processes. In the light of these contributions perhaps 
the problem of flutter analysis as proposed by Kiissner 
warrants further investigation. 

Recently, Goland (reference 3) applied the differential- 
equation method to a uniform cantilever wing and was able 
to carry out the solution of the flutter problem by straight- 
forward methods. In reference 4 Goland and Luke extended 
the solution of the problem of the uniform wing to include a 
uniform wing carrying a fuselage at the semispan and con- 
centrated weights at the tips. Goland and Luke made use 
of the Laplace transform to solve the differential equations by 
operational methods for both the symmetric and antisym- 
metric types of flutter. In both references 3 and 4, the 
objective was to compare flutter speeds and certain flutter 
parameters for specific uniform wings calculated by the 
differential-equation method with the same quantities calcu- 
lated by the Rayleigh method when only the fundamental 
bending and torsion modes were used in the calculations. 
Fairly close agreement between results calculated by the two 
methods were obtained in both references 3 and 4. No 
comparison with experiment, however, was made in either 
case. 
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The results of a systematic series of flutter tests made to 
determine the effect of concentrated weights and concen- 
trated weight positions on the flutter speed of a uniform 
cantilever wing are reported in reference 5. After these 
experiments were finished, the results were compared with a 
theoretical analysis by the Rayleigh method (reference 6). 
In cases where the mass of the weight was of the same order 
as that of the wing and placed so that the distance between 
its center of gravity and. the elastic axis of the wing was a 
considerable fraction of the wing chord, several normal modes 
had to be employed and there was no way of knowing in 
advance just what number should be used. Because of this 
difficulty and because the wing was a uniform wing, the most 
extreme case was chosen from reference 5 and investigated by 
the differential-equation method by following an extended 
procedure of Goland and Luke. The purpose of this report 
is to present the results of this investigation. 

This report consists of the main text and two appendixes. 
In the main text the differential-equation method is set up 
for any uniform cantilever wing with an arbitrarily placed 
concentrated weight and the solution, based on an extension 
of the method used by Goland and Luke, is developed. 
Application is then made to a particular wing-weight system 
used in reference 5, and comparison with experimental 
results is given. The mass of the weight (weight labeled 7a 
in reference 5) was about 92 percent of the mass of the wing 
and at each spanwise weight position the weight was placed 
so that its center of gravity was about 0.41 chord forward of 
the elastic axis of the wing. (It may be mentioned for the 
sake of comparison that, in the numerical example treated in 
reference 4, the mass of the weight was only 39 percent of the 
mass of the wing and was placed 0.1 chord behind the elastic 
axis of the wing.) The geometric aspect ratio of the wing 
was 6, which was considered large enough to warrant the use 
of two-dimensional air forces without aspect-ratio corrections 
for oscillatory instability (not necessarily so for the divergent 
type of instability (see reference 7)). One other simplifica- 
tion was the omission of terms due to structural damping. 
The computed results agree remarkably well with experi- 
mental results, particularly in regard to trends. 

In appendix A the method used by Goland and Luke, 
which includes the derivation of the differential equations, 
for a wing carrying a tip weight is outlined and extended to a 
wing carrying an arbitrarily placed weight. A somewhat 
general but brief discussion of operational methods of solving 
boundary-value problems is included and illustrated with a 
simple example for readers who might be interested but are 
not familiar with the operational approach. 

In appendix B the derivation of the flutter determinant is 
completed, and a method of solving the determinant is 
illustrated by a detailed calculation of the flutter speed for 
the wing and one weight position of the wing-weight combina- 
tion discussed in the text. As a final topic in this appendix 
the solution obtained for the flutter determinant is used with 


the solutions of the differential equations to calculate the 
amplitudes and phase angles of the deflection curves of the 
wing-weight system at flutter speed. 

SYMBOLS 


a nondimensional distance of elastic axis from 

midchord measured in half-chords, positive 
for positions of elastic axis behind midchord 
b wing half-chord, feet 

e x chordwise distance of wing center of gravity 

from elastic axis, positive for center of 
gravity behind elastic axis, feet 
e 2 chordwise distance of weight center of gravity 

from elastic axis, positive for center of 
gravity behind elastic axis, feet 
g gravitational constant, feet per second per 

second 

I mass moment of inertia of uniform wing per 

unit of spanwise length, referred to wing 
elastic axis, pound-second 2 ( mKi ) 

I w mass moment of inertia of weight referred to 

wing elastic axis, foot-pound-second 2 
K ] radius of gyration of wing sections about wing 

elastic axis, feet 

Ko radius of gyration of weight about elastic axis, 

feet 

k reduced-frequency parameter 

L aerodynamic lift force per unit of spanwise 

length 

L y -\~iL y ' =irpb 2 L h 

Le + iL g '=T P b 3 [La-Ln (j +«)] 



l semispan of wing, feet 

location of weight measured from wing root, 
feet 

L h ,L a ,M],,M a aerodynamic coefficients as tabulated in 
reference 8 


M 

aerodynamic moment per unit of spanwise 
length taken about elastic axis 

= 

: irpb 3 

+ 

h|(N 

1 

1 

Me -f -iM$' = 

irpV 

M a — L a ( 2+gT)— 

IF 

weight of wing model, pounds 

m 

mass of wing per unit length 

W, c 

weight of concentrated weight, pounds 

N 

transverse shear force in wing at station x 

T 

torsional moment in wing at station x 

R\, B 2) R 3 

roots of cubic equation 

s 

operator used in Laplace transformation 

t 

time coordinate 


sum of all symmetric polynomial functions in 
R u B 2 , R 3 which are of degree n 
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»0 


V 

V 


bio 

X 

Y(x,t ) 

yM 


EI b 

GJ 


experimental flutter speed for wing without 
weight, feet per second 
flutter speed, feet per second 

reduced flutter speed 

spanwise coordinate measured from wing root 
general mode shape function in bending 
mode shape function in bending after assump- 
tion of harmonic motion (yi(x)+iy 2 (x)) 
flexural rigidity of uniform wing, pound-feet 2 
torsional rigidity of uniform wing, pound-feet 2 


ex — -pij (rn~ f- Ly -j- %Ly ) 

co 2 

P Le-\-iLg') 

CO 2 


K 

p 

A 

Ag 

G(z f t) 

e(x) 

CO 


J 


mass ratio 



air density, slugs per cubic foot 
complex value of determinant 
value of A when real and imaginary parts are 
equal • 

general mode shape function in torsion 
mode shape function in torsion after assump- 
tion of harmonic motion (f? 2 (x) -\-id 3 (x)) 
circular frequency at flutter, radians per 
second 


frequency, cycles per second 



ANALYSIS 


fa) y(O)=y'(O)=0(O)=O • 

(b) EI b y" (l) =EI b y"'(l) = GJd' 0 1 ) =0 

(c) Eh[v”' (!i-0) i+ 0 )] = -^ « 2 [y(h)+e 2 em 

(d) (?J[0 / (Z 1 --O)-0'(Z 1 +O)]=^ tfiea{ii)+K a *o m 

where 


a— -gj- (m +£„+££/) 

P~Eh ^ mei Yl j g-]riLe) 

(mex+My+iMy') 

d=-gj(I+Mo+iMo') 


and where y(x) is the displacement of a chordwise element of 
the elastic axis of the wing at span position x due to bending: 
6(x ) is the corresponding displacement due to torsion; primes 
associated with y and d indicate differentiation with respect 
to x; EI b is the flexural rigidity of the wing; GJ is the tor- 


sional rigidity of the wing; is mass of the weight; m is 


mass per unit length of wing; and io is the circular frequency 
of bending and torsion at flutter. In condition (c) the nota- 
tion y"'(li — 0) indicates that y'"(z) is to have the value 
that it approaches as x^h from the inboard side of the weight 
and y'"(h+0) indicates that y"'(x ) is to have the value that 
it approaches as from the outboard side of the weight. 
Similar meanings are given to 6' {lx — 0) and 0'(Zi + O). 

The quantities Ly+iLy , Le+iLe ' , , and 

Me+iMe' can bo written in terms of tabulated quantities as 
follows: 

L v +iL v ' = irpb 2 L, t 


As mentioned in the introduction the differential equations 
that govern the motion of a uniform wing at flutter speed, as 
derived by Goland in reference 3, and a method of solving 
the equations for a uniform cantilever wing carrying an 
arbitrarily placed weight, based on a method developed by 
Goland and Luke in reference 4, are discussed in appendix 
A. The analysis, therefore, is devoted to a brief discussion 
of the differential equations of motion of the wing, the bound- 
ary conditions, solution of the boundary-value problem by 
means of the Laplace transform, and the solution of the 
flutter determinant. 

The differential equations and boundary conditions that 
govern the motion, at flutter speed, of a cantilever wing of 
length l with a concentrated weight placed l x units along the 
span from the root section and e 2 units forward of the elastic 
axis of the wing, as derived in appendix A, are 


?/ IV (x) — ■ ay {x) ~pe(x)= 0 

(1) 

e" (x) +y y (x) + m(«) = o 

(2) 


Lo-\-iLg' — -rrpb 3 j~L a — 

My-^-iMy' — TTpb 3 — L h 
AI g -f- iMi/ — irpb ' 1 — L a Mn 

Ln (l +a ) ] 

In reference 8 the values of L,„ L a , M,„ and M a are expressed 
in terms of Theodorsen’s F and G functions of reference 9 
and tabulated for various values of the reduced speed v/bio. 

The root conditions (a) and the boundary conditions (b), 
of the boundary- value problem, are the usual conditions that 
must be imposed upon the equations of a vibrating cantilever 
beam (or wing). Conditions (c) and (d) stipulate discon- 
tinuities of determinable magnitudes in transverse shearing 
force and torque, respectively. 
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Applying the Laplace transform (see appendix A) 

e~ sx f(x)dx=J(s) 


to equations (1) and (2) and making use of conditions (a), 
(c), and (d) gives 


and 


s i y(s)-sY 2 -Y 3 +e-* l i[y"'(h-0)-y"'(l 1 + 0)y 
ay (s) — I3d(s) == 0 


(3) 


s 2 d(s) — 01 + e sil [0'(7i— 0) — 0'(Zi+O)] + 50(s) + 77 /( 5 ) =0 (4) 
where 

*W'( 0) 

Y s =y'"( 0) 

6 1 =d f (0) 

Solving equations (3) and (4) for y(s) and 0(s) gives the 
Laplace transform of y(x) and 8(x) , respectively, as 


y(A : 


and 


(s 3 +g5) T 2 + (s 2 +0) Y 3 -\-(j6i — fie — 0)— 0 / (Zi+O)] 


e0) 

(s 2 -\-d)e~ sl '[y"' (^—0) —y'" + + 0)] 
q(s) 


(5) 


QS, s%-ysY 2 -vY 3 -d ia +ye sh [y"'(h — 0) — y"'(h + 0)] . 
^ J 2(s) 


(a~s i )e~° l '[d'(l 1 -0)-d'(l 1 + 0)] 
<l(s) 


(6) 


where 


q(s) =S G + §S 4 — aS 2 -\-y(3 — ab 


Goland and Luke (reference 4) showed that y(x) and d(x) 
could be written as a converging series by expanding the 
transforms (5) and (6) into power series and applying the 
inverse transform. A discussion of this expansion is given 
in section 4 of appendix A where it is shown that 1 jq(s) 
can be written as 


where 


For n7> 3, 


1 1 » T 

To— l 
T^—b 
T 2 =b 2 +a 
Tz= — b 3 — ab — f3y 


5T’ re _i + aF„_ 2 + (ab—@y) T n . 


(7) 


( 8 ) 


When the series expansion of 1 jq(s), equation (7), is sub- 
stituted into equations (5) and (6), the transforms y(s) and 


6(s ) become sums of infinite series with terms of two distinct 
types; that is, terms of types 


and 


A 

S m 

Be~ sx 0 


s 71 


where m is a positive integer. 

A 

The inverse Laplace transform of (see pair no. 3, 

o 

p. 295, of reference 10) for a:>0 is 

r-i M) = Ax” 1 - 1 
\ s m ) (m — 1)! 


(9) 


Be 3X ® 

and the inverse Laplace transform of — (see pair no. 63, 

s 

p. 298, of reference 10) for x^>Xq^. 0 is 

7 -_, $Be~ sx °\ _B(x—x 0 ) m ~ l 
\ s m )“ (m— 1)! 


( 10 ) 


When the expression for 1 fq(s) from equation (7) is sub- 
stituted into equations (5) and (6) and the inverse transforms 
are applied, the following series expressions of y(x) and 6(x) 
can be obtained: 


C co 'T 7 71+2 co y ™2n+4 

co 'J 1 r *% 2 tt-|-3 “1 co T n-f-5 

V 7 M If n a 'Y' — ndj 

(271+3) !_J ' lP £=0 (2w+5)! 


W xA-o)-^ft+o)ig /^;; w . 


-| r rp r 2n+5 


[y"'{U- 0)-y"’(U+m [sg^| 


-h) 2n + 5 


(271 + 5) ! 


- + 


“ T n (x—li) 2n+3 ~\ 
(2ti+3)! J 


(ID 


and 

co / T 7 /y»2n+l co Y 7 ^2n+4 oo '7 7 /y.2«+5 

±nX 


f=o (2w+l)! 


feo (2n+4) ! 


(271+5) ! 


CO 7 1 r 2 n+b 

6ia S (2n+5) 

f +, T n (x-l :) 2 ”+ 5 +, T n {x-k) 2n +n 
L“^o (2n+5) ! n=o (27i+l)! J + 


( 12 ) 


where in both equation (11) and equation (12) the terms 
involving (x—li) are to be considered as zero when x— l x . 

Equations (11) and (12) are general expressions for the 
amplitudes or displacement of a point x of the elastic axis of 
a uniform wing vibrating in bending and torsion under the 
conditions of flutter with an arbitrarily placed concentrated 
weight. When the weight is concentrated at the wing tip, 
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the equations correspond to those obtained by Goland except 
for a difference in root conditions. When the weight is 
concentrated at the root (or if the mass of the weight is 
reduced to zero), the equations reduce to those for a uniform 
cantilever wing. These equations may appear rather formid- 
able in their present form; however, only the first few terms 
of each summation seem necessary for most cases. 

In the derivation of the flutter determinant in appendix B 
it is shown that since terms involving (x— l x ) drop out of 
both equation (11) and equation (12) at x=li, the values of 
y(li) and d(l{) can be obtained from the terms not involving 
(x — . Then, by making use of conditions (c) and (d) 
again, linear expression in Y 2 , Y 3 , and 0j can be substituted 
for the bracketed expressions 

[y'"(h-0)-y'"(h + 0)] 

and 

[0'(h-O)-0'(Z 1 + O)] 

After the substitutions are made, equations (11) and (12) 
will contain only the, three undetermined coefficients Y 2 , Y 3 , 
and 0j for any particular wing-weight system of the type under 
consideration. Observe that conditions (b) have not yet 
been used. If these conditions are now imposed upon the 
equations, there is obtained a system of three linear homo- 
geneous equations in Y 2 , Y 3 , and 0i that may be written for 
reference as 

A t Y 2 +B t Y 3 +CA= 0 (13) 


methods exist, a graphical method was adopted for the pres- 
ent work. For example, a value is assigned to one param- 
eter, preferably vjbw, the flutter determinant is then 
evaluated for this value of v/bu j and several values of the 
other parameter o>. The values of the flutter determinant 
obtained in this manner are complex numbers and if the 
real and imaginary parts of a sufficient number of determi- 
nant values are separately plotted against «, the point or 
points where the real and imaginary parts are equal are 
obtained. If this process for other values of v/bw is re- 
peated, a locus of determinant values with equal real and 
imaginary parts can be plotted against both vjbo) and «. 
When enough points are determined, these plots give the 
values of v/bu and co that cause the determinant to vanish. 

An illustration of the process of solving the flutter determi- 
nant as described in the preceding paragraph is given in 
appendix B, which contains the complete solution of the 
determinant for one weight position of the particular wing- 
weight system described in the section entitled “Applica- 
tion to a Specific Wing-Weight System.” In general, when 
solving the flutter determinant by the preceding method, 
if the assumed values of v/bu and are in the neighborhood 
of their true values, only a few points need be computed to 
obtain a solution. In the absence of experimental values 
of these parameters and in view of the work involved in 
determining other parameters that depend on vjbw, it will 
be found advisable to use simplified methods to obtain 
approximate values with which to start the solution. 


where i= 1, 2, and 3. 

The condition that a system of equations such as equations 
(13) have solutions other than the trivial solution 

Y 2 =Y 3 = 8i = 0 

is that the determinant of the coefficients A i} B t , and C t 
vanish (reference 11). This corresponds to the borderline 
condition between damped (stable) and undamped (unstable) 
oscillations or to the point at which flutter occurs. It will 
be noted that the order of this determinant depends only on 
the order of the system of differential equations. 

The actual coefficients corresponding to A t , B u and C t are 
complex functions of the frequency o>, the reduced flutter 
speed v/bw, and certain determinable characteristics of the 
wing-weight system. The true flutter speed is easily cal- 
culated when corresponding values of a and v/bco are known. 
These quantities may therefore be considered as (the only) 
variable parameters in the determinant of coefficients and 
the problem of finding the true flutter speed is reduced to 
that of finding corresponding values of these parameters 
that cause the determinant, hereinafter called the flutter- 
determinant, to vanish. If v is set equal to zero, the air 
forces drop out and the resulting determinant gives the 
coupled modes of vibration of the wing in still air. On the 
other hand, if « is set equal to zero the nonoscillatory or 
divergence condition is obtained. 

Several ways of solving the flutter determinant are 
mentioned in reference 7. Although more informative 


APPLICATION TO A SPECIFIC WING-WEIGHT SYSTEM 


Attention is now turned to the application of the boundary- 
value problem discussed in the foregoing section to a specific 
problem. The wing-weight system that has been analyzed 
consists of a particular uniform cantilever wing and weight 
combination described in reference 5. The weight was con- 
sidered as concentrated at different specified span positions 
but always at about 0.41 chord forward of the elastic axis 
of the wing. This weight was selected because of its high 
mass compared with that of the wing and because of the 
large eccentricity due to the distance between its center of 
gravity and the elastic axis of the wing. Pertinent data, 
based on measured characteristics of the wing as taken from 
reference 5, with the units in feet and pounds are 


Chord, feet 

Length, feet 

Aspect ratio (geometric) 

Taper ratio 

Airfoil section 

W, pounds 

I, pound-second 2 

EI b , pound-feet 2 

GJ, pound-feet 2 

1/k (standard air, no weight) 
ei, feet 


/3 

4 

6 

1 

NACA 16-010 
3.48 
0.00080 
977.08 
480.56 
32.6 
0.013 


and, based on measured characteristics of the weight, are 


PA, pounds 

e 2 , feet 

I w , foot-pound-second 2 


3.182 
— 0.2728 
0.013625 
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Calculation of the flutter parameters have been made for 
the wing without the weight and for the wing with the 
weight at six different positions. The calculated results are 
compared with experimental results in figure 1 and in table I. 



Figure 1. — Comparison of calculated and experimental flutter speeds for a particular wing- 

weight system. 

TABLE I— EXPERIMENTAL AND CALCULATED RESULTS 
OF FLUTTER INVESTIGATION 


Zi 

(in.) 

C2 

(ft) 

Calculated 

Experimental 

(cps) 

v/boj 

V 

(fps) 

(cps) 

v{b(t) 

V 

(fps) 

0 


25.27 

6.29 

333 

22.1 

7. 22 

334 

11 

-0. 2728 

19.23 

8. 23 

331 

17.4 

8.88 

324 

17 

-. 2728 

28.04 

6.93 

407 

°26.8 

6. 81 

382 

30 

-.2728 

30.68 

8. 18 

526 

(>) 

(M 

(») 

45 

-. 2728 

25.67 

7.45 

401 

O) 

U) 

(«) 

46 

2728 

24.87 

7.06 

368 

21.8 

8. 06 

368 

48 

-. 2728 

23.60 

6. 07 

300 

21.4 

7. 14 

320 


« It is found in reference 5 that good flutter records for this wing-weight system were 
obtained for several spanwise weight positions between the root section and a point 17 inches 
from the root section; but with the weight at 17 inches from the root section the wing appeared 
to diverge. However, the oscillograph records for this case showed two possible flutter 
points: one corresponding to a frequency of 16.3 cps and another corresponding to a frequency 
of 26.8 cps (only the first of these is recorded in reference 5). When the weight was moved 
farther outward from this point, definite divergence was noted until the weight was at a 
point 46 inches from the root section. At this point and from this point to the tip good flutter 
records were obtained. 

1 Divergence. 


- It will be noted in table I that all the calculated flutter 
speeds are within 7 percent of the experimental values and 
the calculated frequencies and reduced speeds are within 15 
percent of the experimental values. The calculated flutter 


speeds are generally slightly higher than the experimental 
values for US17 and slightly lower for ^^46. There is no 
such consistent trend in the other parameters. 

In figure 1 the ratio of both calculated and experimental 
flutter speeds for the wing with a weight to the flutter speed 
of the wing without a weight is plotted against span position 
of the weight. The important thing to note in examining 
figure 1 is that the shape of the theoretical curve follows the 
shape of the experimental curve very closely in the regions 
where experimental flutter was obtained. The horizontal 
dashed line in figure 1 represents the divergence speed for 
the wing as computed by the method of reference 12. Al- 
though the correct divergence speed for different weight posi- 
tions would probably vary, being somewhat lower with the 
weight at the tip than at the root, owing to the effect of the 
presence of the weight on aerodynamic forces, the agreement 
of the approximate value with experimental values is satis- 
factory. 

General expressions for the deflection curves are derived 
in appendix B from which amplitudes and phase angles for 
curves of deflection, slope, moment, and shear in bending 
and amplitudes and phase angles for curves of angular deflec- 
tion and torque in torsion can be computed. The phase 
angles and amplitudes for the deflection and shear curves in 
bending (fig. 2) and the phase angles and amplitudes for the 
angular displacement and torque in torsion (fig. 3) have been 
computed with reference to a unit tip deflection for the weight 
position h = \7 inches. In figure 4 the amplitudes in deflec- 
tion and shear in bending from figure 2 are compared with 
the deflection and shear curves due to the fundamental un- 
coupled bending mode of the wing, and in figure 5 the ampli- 
tudes in angular deflection and torque in torsion from figure 3 



Figure 2.— Plot of amplitude and phase angle of displacement and shear curve in 
bending at flutter for fi=17 inches (amplitude and shear referred to unit amplitude 
at tip in bending). 
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Figure 3.— Plot of amplitude and phase angle of torsional displacement and torque 
for ?i = 17 inches at Butter (amplitude and torque referred to unit amplitude at tip in 
bending). 


are compared with the angular deflection and torque curves 
due to the fundamental uncoupled mode in torsion. There 
is a notable difference in the shape of the amplitude curves 
computed by the present method and those computed from 
the first normal modes. This discrepancy indicates that 
several modes would have to be employed to obtain satis- 
factory results by the Rayleigh type analysis. 

CONCLUDING REMARKS 

The method discussed in this report is not limited to a uni- 
form cantilever wing with a single weight. By proper atten- 
tion to the boundary conditions the theory can quite easily 
be extended to apply to a uniform wing carrying any number 
of arbitrarily placed weights, one of which might be consid- 
ered as a fuselage and made to yield the so-called symmetric 
and antisymmetric types of flutter. Furthermore, for con- 
venience of application, theoretical values of two-dimensional 
air forces have been used. However, since the method does 
not depend on the particular form of air forces involved, any 
known or available aerodynamic data could be used. In 
any event, the method is tedious and would, therefore, 
not be recommended over the Rayleigh type analysis when 
it might be known that only the first few normal modes of 
the structure are sufficient to give satisfactory results. 

For wings that are not uniform the differential equations 
for flutter conditions reduce to ordinary differential equations 
with variable coefficients. In this case the solution would, 
in general, be much more difficult to obtain. For general 
cases there would be no advantage in the operational method 
of solution although an iterative process probably might be 
used to great advantage. 

In conclusion it is pointed out that the location of such 
masses as engines, landing gears, and fuel tanks might be 
used to advantage in increasing the flutter speed of a given 
wing* s h° wn by the particular problem analyzed herein 



Figure 4 .— Plots of amplitudes in bending displacement and torque and the corresponding 
curves computed for the first uncoupled normal mode in bending for fa = 17 inches (ampli- 
tude and shear referred to unit amplitude at the tip in bending). 



Figure 5. — Plots of amplitudes in torsional displacement and torque and the corresponding 
curves computed for the first uncoupled normal mode in torsion for Zr =17 inches (ampli- 
tude and torque referred to unit amplitude at tip in bending). 

and by other experiences a definite region exists, peculiar to 
a given wing, in which masses added forward of the elastic 
axis of the wing tend to increase the flutter speed of the wing. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November SO, 1948. 







APPENDIX A 

OUTLINE AND EXTENSION OF METHODS OF FLUTTER ANALYSIS AS PRESENTED IN REFERENCES 3 AND 4 

and at ( x-\-dx ) the nose-up torsional moment is 


1. DERIVATION OF THE DIFFERENTIAL EQUATIONS THAT GOVERN THE 
MOTION OF A WING AT FLUTTER SPEED 

Consider a spanwise element of incremental length dx at 
station a; of a wing oscillating in bending and torsion in a free 
stream of fluid (see sketch 1). 



Station x 


The displacements F and 8 of an element of the elastic axis 
are functions of x and t. In order that this clement remain 
in dynamic equilibrium the external forces and moments on 
the element must balance the inertia forces and moments. 

The external forces and moments consist of transverse 
shearing forces and torsional moments, which are transmitted 
from one element of the wing to the next, plus the aero- 
dynamic lift force and pitching moment and internal or 
structural damping. Structural damping is not taken into 
consideration in this discussion, although its inclusion would 
add no computational difficulties. 

The transverse shearing force acting upward at x is 


N=-EI, 


d 3 F 

d 3 x 


(Al) 


and that acting downward at ( x+dx ) is 

dA r , Tn t d 3 F 

da; 


N+ dx— —El b ^ —EI b hh dx 


d 4 F 

da; 4 


(A2) 


Similarly the nose-down torsional moment acting at x is 

d0 


T— GJ 


dx 

d 2 F 


(A3) 


d 2 e 




d 4 F 


T+ dY dx=GJ dx + GJ dY dx 


(A4) 


The two-dimensional aerodynamic forces acting on an ele- 
ment dx of an oscillating airfoil have been derived by 
Theodorsen (reference 9) and can be written as a lift force 
and aerodynamic moment acting about the elastic axis of 
the wing, respectively, as 


L dx—^oi 2 L y Y-\-uL y ' 


dx (A5) 


Mdx=( orM y Y+wMy' ^ + <YMqQ -}- uMg' J 


dx (A6) 


The inertia force of the element dx can be written 


' d 2 F . d 2 0 N 

m ~w +mei W y 


and the inertia moment as 

d 2 8 


r ^ ^ d 2 F\ 

I -W +me ‘ w) 


dx 


dx 


(A7) 


(AS) 


Diagrams of the forces and moments acting on an element 
of wing of length dx at station x are shown in sketch 2. 




Imposing the conditions of dynamic equilibrium of the 
element at x by equating inertia forces to external forces and 
inertia moments to external moments gives the two dif- 
ferential equations that govern the motion of the wing: 


/ dF 


-EI b ~+a 2 L y Y+aL v ' d --+ a> 2 £,e+«2* 


, de 
dt 


I W +me ' W =ej 5J + “’ M ’' 7+ “ Mt ' d it 


(A9) 
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2. BOUNDARY CONDITIONS FOR A UNIFORM CANTILEVER V.ING CARRYING 
AN ARBITRARILY PLACED WEIGHT AT FLUTTER SPEED 

The boundary conditions that must he imposed upon 
equations (A9) for a uniform cantilever wing are 

(1) F(0,i)=0 

( 2) , 

(3) ©(0,2) — 0 

(4) Eh\j^Y{x,t)^ x _ =0 

(5) EI b ^^Y(x,t)~^ x= =0 

( 6 ) Gj[le(x,t)] x = 0 

These are the usual conditions that must be imposed on a 
vibrating cantilever beam. Condition (1) is the condition 
that the end at x=0 is supported (either hinged or built in). 
Conditions (2) and (3) imply that this end is fixed or built 
in. Conditions (4), (5), and (6) imply, respectively, that 
there is no bending moment, transverse shearing force, or 
torsional moment acting at the tip x—l. 

If there is an arbitrarily placed weight on the wing, other 
conditions must be imposed that will determine the effect 
of the weight upon the motion of the wing. If the weight 
is considered as concentrated at some point on the chord 
line at station x=li, it will create discontinuities in both 
transverse shear and torsional moment. The magnitude of 
these discontinuities are known functions of the mass of the 
weight, the location of the weight, and the acceleration of 
the wing. The remaining conditions required to complete 
the boundary-value problem for the general motion of the 
weighted wing are, therefore, 


the frequencies in bending and torsion are equal. Therefore, 
only the particular form that the solution to the boundary- 
value problem has when these conditions prevail need be 
sought. These conditions imply that Y(x,t) and Q(x,t) are 
of the forms 


Y(x,t)=y(x)e ia>t \ 
Q(x,t) — 9(x)e iat ) 


(A10) 


where, on the right-hand side of equations (A10), y and 6 are 
now complex amplitude functions of the span coordinate x 
from which the shape and phase relation of the wing at 
any fixed time during flutter can be obtained. 

If the values of Y and 6 from equations (A 10) are sub- 
stituted into both differential equations (A9) and into the 
boundary conditions, the problem is greatly simplified. The 
differential equations become independent of t and appear 
as ordinary differential equations with constant coefficients. 
After making the substitution and rearranging terms, the 
equations of motion can be written as 


EI b L y -\-iLy' ) <s?y — (mei-j- Le J riLg')<x> 2 d=0 

d 2 R 

GJ ^ + (me 1 +M J ,+fM/)w 2 y+(/+M 9 +tAf 8 ')co 2 6i=0 


Y (All) 


or more simply as 


dSj 

dY 


ay — /30=O 1 


(Pe 

dx 2 


~\~yy ~r 50= 0 




j 


(A12) 


The boundary conditions also become independent of t and 
can be written as follows: 


(!') 2 /( 0 ) =0 


(7) EI t \ fS Y(x, i)l -fa Y(x, ()1 

( L ox Ji=(! r o) L ux Jt=(i,+o) 

= ~hr [* sp *>+*>' w 6 <*' «]„ ,, 

For the purpose of flutter analysis it is assumed that the 
motions in both bending and torsion are harmonic and that 


(2') y'(0)=0 
(3') 6(0) =0 
(4'j y"(l)=0 
(5') y"'(l)= 0 
(6') 0'(O= 0 

(7') EI b [y'''(l l -0)-y''\l 1 +0)] = -^a, 2 [y(l>)+e 2 d(l l )] 
(S') GJ[6'(U-0)-d'(U+0)}^-- w - JfaydJ+KMh)] 
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3. SOLUTION OF BOUNDARY- VALUE PROBLEMS IN ORDINARY DIFFEREN- 
TIAL EQUATIONS BY OPERATIONAL METHODS AND APPLICATION TO A 
BEAM CARRYING AN ARBITRARILY PLACED WEIGHT 

The boundary-value problem given by equations (A12) 
and conditions (1') to (8') can be solved by straightforward 
methods of solving ordinary differential equations with 
constant coefficients. The operational method, however, 
is a much easier and shorter approach, particularly in view 
of the discontinuities in shear and torque. 

Briefly, the solution of a boundary-value problem by 
operational methods consists of applying the Laplace trans- 
form to the differential equations, the initial conditions (root 
conditions when applied to beam problems), and certain 
forms of other boundary conditions; of solving the resulting 
system for the transform of each dependent variable; and 
then applying the inversion integral to the results. The 
remaining boundary conditions are then used to set up 
relations among whatever undetermined parameters that 
might remain. 

In the case of flutter analysis a complete solution to the 
equations is not needed but only the conditions under which 
an unstable equilibrium may exist. The relations that can 
be set up between the undetermined parameters correspond 
precisely to this condition In other words these relations 
appear as a system of homogeneous equations and the satis- 
faction of the condition that this system of equations have a 
common solution other than the trivial solution corresponds 
to the borderline condition separating the damped and 
undamped oscillations of the wing. 

The Laplace transform of j{x) is 

L{f(x)}*= ^ e~ sx j(x) dx=J(s) (A 13) 

where s may be real or complex and z)>0. The sufficient 
conditions that this infinite integral exist are that f(x) have 
no infinite discontinuities for and that f{x) be of ex- 
ponential order as °° . (See reference 10.) In other words 
finite discontinuities such as those appearing in the foregoing 
problem do not invalidate the operational approach. 

The Laplace transform of the nth derivative of a con- 
tinuous function with continuous derivatives, for which the 
function and all its derivatives are of exponential order, 
can be obtained directly from equation (A 13) as 

L{f n (x ) } =s n f (s) — s” -J /(0) — s” -2 /'(0) — . . . — /" -1 (' 0) 

(A 14) 

The Laplace transform is linear in the same sense as 
differentiation or integration. That is, if a, and b { are 
constants 

Z{a n / n (x)+a B _i/” _1 (z) + . . . +a o j(x) + b m 0 m (x)i- 
. . . -\-b o 0(x ) } —a n L{j n {x ) } J ra n _iL{f n ~ l (x ) } -f- ...-f- 
a 0 L{f (x)}+b m L{d m (x)}+ . . . +b 0 L{d(x)} (A15) | 


Thus the Laplace transform of a linear differential equation 
with constant coefficients is generally a sum of expressions 
similar to equation (A14). 

In equation (A14) the quantities /(0), /'(0), . . . , 

/” -1 (0) are the boundary conditions at the origin of the 
dependent variable (wing root) that corresponds to constants 
of integration. When these quantities are given, they are 
put directly into the transformed equation. When the 
quantities are not given, they correspond to what has been 
called undetermined parameters in the preceding paragraphs 
and must later be determined in terms of other boundary 
conditions. 

Finite discontinuities in a function or any of its deriva- 
tives are taken into account by proper attention to the 
limiting values that the function or its derivatives have on 
the two sides of the discontinuity. In particular, if a func- 
tion and its first n derivatives are of exponential order, if 
the first (r> — 2) derivatives are continuous, if the (n— l) sl 
derivative has a finite discontinuity at x 0 , and if the nth 
derivative is continuous except for a singular point at x 0 , 
(see sketch 3), the Laplace transform of the nth derivative 
has the form 

L{f n (x)}=sJ(s) — s n ~ l f(0)— . . . — 

s/"~ 2 (0) — / n_1 (0) — e _s:co [/ n_1 (ro+0) — ; J n ~ l {x Q ~-0)\ 

(A16) 

where /(x 0 +0) is the value of/(x) as x approaches x {) from the 
right and f(x Q — 0) is the value of f(x) as x approaches x Q 
from the left. In other words the terms in the brackets 
express the magnitude of the discontinuity in/" _I (x) at x 0 in 
the (n— 1) SJ derivative at x 0 . 



An examination of the boundary-value problem, equation 
(A12), shows that the transform will be given by a sum of 
expressions precisely of the form of equation (A16). 

In order to interpret the transformed function f(s) in 
terms of the original function /(x) , use may be made of the 
inversion integral discussed in text books on operational 
calculus; or one may refer directly to tables of transform. 
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As a simple example the operational method is applied to 
a cantilever beam carrying an arbitrarily placed weight and 
assumed to be vibrating in a vacuum in bending only. 

The boundary-value problem for this case can be written 

EI t (Al7) 

(a) y( 0 )= 2 /'( 0)=0 

(b) .y"(l)=y'"(l) = 0 I (Alg) 

(c) EI,[y'"(h-0)-y"'(.l,+0))—^ y(l t ) 

y j 

where the symbols have the same meaning as in equation 
(A12). 

If the root conditions (a) and the boundary condition (c) 
are used, the transformed problems solved for y(s) gives 


computed and are plotted in figure 6. Corresponding re- 
sults have been computed by a 20-station process of iteration 
discussed in reference 13 and plotted in the same figure. 



Figure 6.— -Plots of deflection and shear curves computed from the first uncoupled modes 
in bending by the differential-equation method and by the 20-station iteration process of 
reference 13 (referred to unit tip deflection) . 


i. REPRESENTATION OF THE INVERSE TRANSFORM OF THE BOUNDARY- 
VALUE PROBLEM, EQUATION (AI2), BY A POWER SERIES 


?/ 0 ) 


sY 2 , Y z 

S i — a 4 _ 1 _ S 4 — a 4 


W w u> 2 y(h) s i x 
gEI b s 4 -« 4 


(A 19)' 


where, for brevity, Y 2 =y"( 0), Y z =y"' (0), auda 4 =^' 


The transform of both y(x ) and d(x) of equation (Al2) is 
of the form 


/(s) 


PM | P 2 (s) 

M ^ q(s) 


e~* o 5 


(A22) 


The inverse transform of equation (A19) is (see pair nos. 
31 and 32, p. 296, and relation 12, p. 294, of reference 10) 

Y Y 

y(x )—~2 (cosh ax — cos ax) -f 2 ~s (sinh ax — sin ax) + 


where P\{s) and P 2 (s) are polynomials both of lower degree 
than q{s). Neither P x (s) or P 2 (s ) have common factors with 
q(s) where q(s) is of the specific form 

q(s) =s 6 +as 4 -J-&s 2 +c= (s 2 — R x ) ( s 2 —R 2 ) ( s 2 —R 3 ) (A23) 


W w u 2 

2a 3 gEI b 


?/(Zi)[sinh a(x— 1\) — sin a[x—l x )] 


(A20) 


or 


y( x )~ 2 ^ ( cos h aX — cos aX ) 'T (sinh ax — sin ax) -j- 

W w a , 2 rY 2 , x. , 

2 a 3 gEI b [2a 2 ( ° 0sh a/l COS + 

F 3 . “1 

2^5 (sinh all — sin al x ) [sinh a{x— l x )~ sin a{x—l x )] 

(A21) 


where the last bracket is zero when x—h^Q. 

Imposing boundary conditions (b) gives two homogeneous 
equations in Y 2 and Y 3 . Each value of a that will cause the 
determinant of the coefficients of Y 2 and Y 3 to vanish corre- 
sponds to a mode of vibration. 

This result has been applied to the wing and weight dis- 
cussed in the text of this report with the weight located 17 
inches from the root. The deflection and shear curves due 
to the first uncoupled modes in bending only have been 


where the coefficients a, b, and c and the roots squared B x , 
R 2 , and R 3 are complex. The inverse function associated 
with such a transform gives f(x) in terms of circular and 
hyperbolic functions of x^Ri, but with the results in this 
form the process of solving the flutter determinant becomes 
very cumbersome. 

By making use of the properties of symmetric functions, 
Goland and Luke (reference 4) outlined a simple method of 
obtaining the multinomial expansions or Maclaurin’s series 
for the transforms of equations (Al2). The inversions of 
these expansions give y{x) and d(x) in the form of convergent 
series and thus circumvents the meticulous task of finding 
the roots of q(s). 

For the development of these series it is first necessary to 
consider q(s ) as a cubic in s 2 ; namely, 

ff(«) = n (s a -fi t )=-« 6 ^(l-p*) (A24) 

By making use of the binomial theorem, 1 /q(s) can be 
written as 


S !( l)“?,?,( 1 + ? + V " + l ® + ’ ’ <A25) 
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Equation (A25) is independent of any interchange of the 
parameters R u R 2 , and R 3 and thus satisfies the description 
of a symmetric function in these parameters. (For a dis- 
cussion of symmetric functions see reference 11 or any text 
on higher algebra or theory of equations.) If the indicated 
multiplication in equation (A25) is carried out, the results 
can be written 

^ = ^(r 0 +-p+^ 2 + . . . +^+ . . .) (A26) 

where the general term T n represents the sum of all possible 
symmetric polynomials in R u R 2 , and R 3 which are of degree 
n and with all coefficients unity. By making use of Newton's 
identity relative to symmetric polynomials, that is 

T n — — aT n _i — bT n _ 2 — cT n _ s (A27) 

where the value of any Tn-j is to be disregarded when 
n— y<C0, every T n can be written in terms of the coefficients 


a, b , and c of equation (A23) ; for example, 
T 0 =l 
jTi= —a 
T 2 =a 2 ~~b 
T 3 =z—a 3 +2ab—c 


(A28) 


With the aid of equation (A26) and equations (9) and (10) 
of the text, the inverse transform of equation (A22) or of 
y(s) and 6(s) can therefore be written as a sum of terms of the 
type given in equations (9) and (10) where the T n ’s enter as 
coefficients in the numerator and are easily evaluated in 
terms of the coefficients of a known cubic equation. In the 
application to flutter analysis only the first few T n ’s are 
usually necessary because the resulting series is generally 
found to be highly convergent. 


APPENDIX B 

DERIVATION OF THE FLUTTER DETERMINANT AND SAMPLE CALCULATIONS 

INTRODUCTION 

In this section the flutter determinant is formally derived and the method described in the text for solving the determi- 
nant is illustrated with sample calculations for a specific example. Also final expressions for the deflection curves are given 
from which amplitude and phase-angle curves of deflection, shear, and torque are calculated for a specific case. The calculated 
amplitudes are compared with corresponding curves computed from the fundamental uncoupled modes in bending and torsion. 

DERIVATION OF THE FLUTTER DETERMINANT 

In equations (11) and (12) of the text it is first necessary to evaluate the expressions 

[2/'"(h-0)— ?'"&+<))] 

and 

ffi (h — 0) — 0' {li -{- 0) ] 

in terms of Y 2 , Y 3 , and 0 U Since terms involving (x—k) drop out of both equation (11) and equation (12) for x—lx , the values 
of y(h) and 6(h) can be obtained directly from these equations. The values of y(k) and 0(Zi) substituted into conditions (c) 
and (d) of the text give the desired relations; namely, 


and 


y , "(h-0)-y"'(l 1 +0) = [y{U)+e 2 e(h)} 


EI b g 


- » rp ] 2 n- f 4 co T / 2 n+2 ~\ 


+ 


1 + 


r o, r r ] 2n+5 o: T 1 2 n+i~\\ 

d 1 ifi-eta) S 4~~~r y + e 2 
L nTo (272.-1-5) ! t=0 (2/i -f- 1) !J) 


ww 


e' ( l, - 0) - d' (l, + 0) - ^ [e 2 y (/,) + KMh) ] 


g 1^45!+^ g (i^yy]+ 
n [to-Bfy) g ^ +. c g ^5,]+ 

r co 7 1 7 2n+5 . T J, 2n + l ~\) 

* [to-KM g §^y-!+^ S ^pryi] 


(Bl) 


(B2) 



where 


FLUTTER OF A UNIFORM WING WITH AN ARBITRARILY PLACED MASS 
Substituting equations (Bl) and (B2) into equations (11) and (12) gives 

y(x) *=ht(z) Y 2 +h 2 (x) Y s +h 3 (x)6 l 
6 (*) =(J\ (x) Y 2 +g 2 (x) Y 3 + g 3 (x) 0 X 
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(S3) 

(B4) 


^ T nX + wwr T n i x 2n+i ^ w^-ir » T n (x-hy n + 5 

{ g (2n+2)! + n=Q {2n+Y)\ + Ehg [_ (5 27) g (2r7+4)! + g (2 t 7+2)!J [_ 5 S (2?7+5 


T n (x~l\) 2n+ * 


(2*7+5)! (2*7+3) 


n+3-l 

rj~ 


m^ 2 

GJg 


[ 


(e 2 S-K 2 2 ) S 


+ T n h 2n+i 


=o (2n+4)!" rg2 S ( 2 n+ 2 )!jj^ (2*7+5)! 


W" +2 1 ^ T n (x~hy n +> 


LW=a V- j y n+5 4-^ y» g 2n+ 3 .a W n+5 L ^ ^i 2n+3 1 T n {x-hy n+5 , ^Cr-ZO^H 

U g (2*7+5) ! + g (2*7+3) ! + #/ s 0 [_ 5 27) S (2*7+5)! + S (2*7+3) J |_ 5 S (2*7+5)! + g (2*7+3)! 


w ” +5 I ^ W* +3 1 ^ T n {x-hy n 
GJg L ( 2 2T) S(2n+5)! +e2 S(27i+3)!jS (2?i+5)! 


/ \2n+S 


r ^ 2ra+5 | w «,Q> 2 r ( , ^ w n+6 . . v a w. 2 ” +l ir^ 

2^ (2?l+5)! + ^ |_(/3— e 2 «) g (27l+5)! +«2 g (2^Fl)lJ [ 5 g (271+5)1 


~^) 2 ”+ 5 , ^ T n {x-k) 


n=o (277+3) 


2n+3-j 

jr“J“ 


/3PF, 


gw [w-*™ a (oti+^ ! a oti] a ^ 


(x-Z,) 2 ^ 5 


(271+5)! 


= — v V* ^+_%T^_yu ^ , +\ T&"*' “lr +, T n (x-ky n +* ^ T„(*-«**+H 

PU 7 g (2*7+4) ! + &/<? [_ ( - 2 7) g (2*7 + 4) ! + * 2 g (2*7 + 2) !_![_' S (2 t 7 + 5)! g (277+1) ! J 


~ET^ [_ ( * 27) g(277+^' + S - 


(277+4)1^^0 (277+ 2)!. 


s 


Tjx-uy^ 


n=o (277 + 5)! 


g 2 (x) = ~y^, 


T n x 2n+S , W w « 2 


(277+5)! GJg 


[ 


® f 1 2 n-f 5 rp i 2n +3 

(e 2 8-K 2 2 y) 23 S " 1 


w [ (5 - e27) a mi+g Mi] a 


^0 (277+5)! 1 02 ^ (277 + 3) 

T„u 2n+z "] » T n (x-uy n+h 


ii r«s 

U L n= 0 


Tjx-uy *+ 5 - r+x-z,) 2 »+H 


(277+5)! 


(277+1)! J' 


(277 + 5)! 


and 


„ 7 1 r 2n + l m V> ™2n+5 TT/ 

j3(*) = S7^nvr«S ” ' 1 


if+o (277+1)! ^ (277 + 5) 


w + r ® T Z, 2B + S 

i+wL^-^S&yi 


^2 2 S 


r»Zi 2n+1 

i+o (277+ 1) ! 


J L «=o v 


(as— W 


2n+5 


-s 


(277+5)’ 

T„(x— Z I ) 2B+1 * 


2it+5 


> n+1— 1 

(277+5)! (277+1)! J“ 

^- 7 f (R e,«) y" r ^ 2n+5 -|. c V 11 ^ n(x-ZQ 2 ” 

EI b g 2 ^ (2/7+5) !^ 02 g (277+1) !jg (2-77+5)! 

By imposing conditions (b) of the text 

y' , (i)=y'"a)=e'(l)= 0 

upon equations (B3) and (B4), three equations are obtained (written in the text as equation (13)): 

A.iY 2 -\- B iY 3 -\- C i 6 x =Q 

where i— 1, 2, and 3 and 

Ai=hi'{l) C x = h 3 "(l) 

A 2 =W"(l) B 2 =h 2 "'(l ) C 2 =h 3 "'(l ) 

A 3 =g/(l) B 3 =g 2 '(l ) C 3 =g 3 '(l) 

Imposing the condition that the equations (13) have a solution other than the trivial solution Y 2 =Y 3 —0 1=0 results 
in the flutter determinant 


A — 


A x 

B x 

c\ 



A 2 

b 2 

c 2 

= 0 

(B5) 

A 3 

B 3 

Os 
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SAMPLE CALCULATION OF FLUTTER SPEED AND DEFLECTION CURVES 

A method of solving the flutter determinant given in the text is illustrated here by the solution of the determinant for 
the wing-weight combination discussed in the text when the spanwise location of the weight is 17 inches from the root. The 

V \ 

values of that are chosen are in the neighborhood of the experimental value and have available tabulated values of 

Theodorsen’s function C(k) —F-\-iG. 

V 

Table II shows the actual computations required to evaluate the coefficients A u B i} and C t for ^=7.1429 (Jt=0. 14) 
and two values of ^ =/ (f—25 cps and/=28 cps). From colums @, (38), and © the determinant for f—25 cps is 


or 


A— 


(14.9200— 2. 8574f) 
(11.8000— 3.6695i) 
(0.17030 — 0.66 134f) 


(12.8320— 2.031 5f) 
(10.2970— 2. 8566i) 
— (0.09077+0.59341f) 


— (7.3286— 0.6002K) 

— (5.4711 -0.93233f) 

— (0.41138— 0.28864i) 


Similarly, for/— 28 cps, 


A= 


or 


(18.6380 — 3.8115i) 
(15.5930 — 5. 0935f) 
(0. 04177 + 0. 87098i) 


A=1.0326 - 0.6948i 

(15.0860 — 2. 6399i) 
(13.0080-3.7946f) 
— (0.23526+0.759480 

A= —0.4029 - 0.0312i 


-(9.1238-0.854330 

— (7.1158 — 1.39880 

— (0.51403— 0.370170 


The determinant was evaluated in this manner for the same value of v/b<o and several other values of/. The process was 

then repeated for j- = 6.25 and several values off and for j— = 5.00 and several values of/. The real and imaginary parts of 
* Deo Oco 

the evaluated determinant for each value of vjbco and the corresponding values of / are separately plotted in figure 7. The 
ordinates of the intersections of the different pairs of curves of real and imaginary parts were scaled in figure 7 and plotted as 

A e against both v/bco and / in figure 8. The zero ordinates of these curves give the value of v/bco (j^=6.93^ and the values of 

/ (/= 28.04 cps) for which the determinant vanishes. From these values the flutter speed is readily calculated to be 


v— (bco) (6.93) = (2t rbf) (6.93) 


(2+) (28.04) (6.93) 
3 


407 fps 


As pointed out in appendix A the deflection curves at any specified time are given by equations (A10) 

Y(x,t) =y(x)e iut —y(x) (cos wt+i sin cot) 

G(x,t)=6(x)e iat = 6(x) (cos co£+f sin cot) 




(a) A. against frequency. (b) A, against reduced speed. 

Figure 8. — Plots of A, against frequency and reduced speed, n =17 inches, 


Figure 7. — Plots of A against frequency for particular values of reduced speed, =17 inches. 
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where final forms of y(x) and 6(x) are given by equations (B3) and (B4) and where, at least, the relative values of the unde- 
termined coefficients Y 2 , Y 3 , and 0 y in equations (B3) and (B4) must be known. If the set of values of vfbco and o> that satisfy 
the flutter determinant is used to determine the coefficients A { , B it and C y in equations (13), there is obtained a system of 
three homogeneous equations in the three unknowns Y 2 , Y 3 , and 6 y that have solutions other than the trivial solutions 
Y 2 : =Y 3 =0i=O. If these equations 'are each divided through by any one of the unknowns, say Y 2 , there is obtained a consist- 
ent system of three equations in the two ratios YJY 2 and 6JY 2 . Any two of the three equations can therefore be solved for 
these ratios. Consequently, equations (B3) and (B4) can be written with one undetermined parameter that appears as a 
factor in each equation. Furthermore, since the coefficients A t , Bi, and Ci are complex numbers, the ratios Y y /Y 2 and 6 y jY 2 
are complex numbers and equations (B3) and (B4) contain complex coefficients. The real and imaginary parts of these 
equations can be separated and the equations written as follows: 


y(x) = Y 2 [y 1 (x)+iy 2 (x)U 
0(x) = Y 2 [e 2 (x) +ie 3 (x)] ) 

If these relations are substituted into equations (AlO), 

Y(x,t) ~Y 2 {yi(x) cos cot—y 2 (x) sin w£+t[y 2 (z) cos cot+y^x) sin cot ] } 

0(x,t) = Y 2 {6 2 (x) cos ut—d 3 (x) sin wt-\-i[d 3 (x) cos coi+0 2 (^)sin w£]} 
or 


(B6) 


(B7) 


where 

and 


Y{x,t) = Y 2 -yj[y y {x)f+[y 2 {x)] 2 [cos (« i+^O+i sin («*-}- <fo)] 
9(^,0 — Y 2 V[^(a:)] 2 +[0 3 (a)] 2 [cos M+02)+4 sin (& >t+<f> 2 )] 


(BS) 


0i = tan 1 


y»(g) 

yiG*D 


0 2 =tan~‘ 


d 3 (x) 
02 (*) 


and where <j>i — <j> 2 represents the difference in phase angle between bending motion and torsion motions at x. 

The real parts of equations (B8) are interpreted to mean the motions in bending and torsion taken in a positive sense. 
The imaginary parts can then be interpreted as representing these same motions with a phase shift of 7r/2 radians. 
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TABLE II— SAMPLE CALCULATIONS OF COEFFICIENTS A { , B { , AND C< 

[fc=0.14] 


® 

© 

© 

© 

© 

© 

(cps) 

O' 

P 

y 

S 

^ W n+2 

o (2n+2)! 

25 

0.65247 — 0.20341* 

-0.48470 + 0.02865* 

0.03178+0.05169* 

0.16724-0.03661* 

0.98302+0.00375* 

28 

0.81846 — 0.25515* 

-0.60800 + 0.03594* 

0.03986+0.06485* 

0.20984-0.04593* 

0.97789 + 0.00468* 


© 

© 

© 

® 



^ rjr 2n+4 

(2r* + 4)! 

[©X © + ©]F 2 

^ TJ, 2 " +S 
(2** + 5) ! 

^ Tnil 2 " +3 
^ (2** + 3) ! 

[©X® + @]F 3 


0.16622+0.00033* 

(1.0108 — 0.00228*) F 2 

0.04721 + 0.00007* 

0.46744+0.00127* 

(0.47534-0.00046*) F 3 


0.16582+0.00041* 

(1.0127 — 0.00285*) F 2 

0.04713 + 0.00009* 

0.46583 + 0.00157* 

(0.46682-0.00058*) F 3 


© 

© 

® 

© 



[® X®]0i 

y(h) =©+©+© 

[~©X©]F 2 

[-®X®]F 3 

[— ©X®]0i 


( — 0.02289 + 0.001 32*)0i 

(1.0108-0.00228*) F, 
+ (0.47534-0.00046*) F 3 
+ ( — 0.02289 + 0. 00132*)Si 

(-0.00527-0.00860*) Y 2 

(-0.00149-0.00244*) F 3 

(- 0.03082+ 0.00956*) 0, 


(- 0.02866 + O.OO164*)0, 

(1.0127-0.00285*) Y 2 
+ (0.46682-0.00058*) F 3 
+ (— O.O2869+O.OO164*)0i 

(-0.00658-0.01077*) F 2 

(-0.00187-0.00306*) F 3 

( — 0.03859 + 0.01 195*)0i 


© 

@ 

© 

© 

© - 


^ M"+‘ 
(2ra + l)! 

[® + ©]0i 

6(1 1) =® + © + @ 

@+e 2 ® 

W CO 2 

7B7® 


1.3695 + 0.00727* 

(1.3387 + 0.01683*)^ 

(-0.00527-0.00860*) F 2 
+ (-0.00149-0.00244*) F 3 
+ (1.3387 + O.O1683*)0, 

(1.0122+0.00007*) Y 2 
+ (0.47575+0.00021*) F 3 
+ (- 0.38791 -O.OO327*)0i 

(-2.5332-0.00017*) F 2 
+ (- 1.1907-0.00051*) Fj 
+ (O.97O82+O.OO818*)0i 


1.3579+0.00899* 

(1.3193 + 0.02095*) 0i 

(-0.00658-0.01077*) F 2 
+ (-0.00187-0.00306*) F 3 
+ (1.3193 + 0. 02095*) 0i 

(1.0145 + 0.00008*) Y 2 
+ (0.46733 + 0.00026*) F 3 
+ ( - 0.38839 - 0.00407*) 0i 

(-3.1848-0.00026*) F 2 
+ (- 1.4671-0.00081*) F 3 
+ (1. 2193 + 0. 01278*’)0i 
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@ 

© 

© 

© 

© 

02© + Ki*® 

W w 2 _ 

y 57 ® 

^ FJ 2 " 

*-• 2w! 

_ T 72 71+2 

^ (2n +2)! 

^ TJ*-« 
f-> (2n+3)! 

(—0.27634— 0.00056i) Y 2 
+ ( — 0.12982 — 0.00021i) F 3 
+ (0.19038+ O.OO195i)0, 

(■ -1.4062- 0.00286i) Y. 

+ (-0.66059-0.00107i) V 3 
+ (O.96875 + O.OO995i)0i 

7.0385-1.9729i 

10.0478 — 0.77851i 

11.424-0.356541 

(- 0.27704- 0.00070i) Y 2 
+ ( - 0. 12755 - 0.00026t) F 3 
+ (0.18928+0.002431)0! 

(- 1.7683 - 0.00449i) Y 2 
+ (-0.81415-0.001680 Y 3 
+ (1.2O82+O.O1553i)0, 

8.7043-2.5064i 

10.6050-0.97996i 

11. 636 -0.44942i 

© 

© 

@ 

© 

© 

^ TJ 2 »+i 
*-> (2» + l)! 

^ T n (l-k) 2 ^ 
. ^ (2n + 3)! 

©X© 

©X© 

^ T n (l—l 1 ) 2n+l 
~ (2n+l)! 

8.0051 -1.3737* 

2.8156 + 0.003211 

-1.3648+0.07913i 

0.47114-0.10255i 

2.7457-0.09258i 

9.0919-1.73201 

2.8018+ 0.00375; 

-1.7036+0.09843i' 

0.58810-0.127891 

2.7925-0.116921 

© 

© 

© 

© 


[@+©X@] Yi 

+ [© + ©] ^’3 + [®]^l 

-[®X@X@1 

+1©X@X@] 

+[®X®]=»"(0 

^ (2n— 1)! 

^ T n (l—h) 2n + 2 
^ (2n+2)! 

®X© 

®X® 

(14.9200-2.8574i) Y 2 
+ (12.8320— 2.0315i) Y 3 
+ (-7.3286 + 0.60021 i)0, 

7.0966-2.32341 

3.3038- 0.01793i 

-1.6008 + 0.103361 

0.55204-0.123961 

(18.6380- 3.81 15i)F 2 
+ (15.0860-2.6399i) Y 3 
+ ( — 9.1238+ O.85433i)0i 

9.1581-3.01441 

3.2979-0.02306i 

-2.0043+0.13256i 

0.69097-0.15630i 

© 

© 

@ 

© 

® 

^ T„(l—k) 2n 
2 n! 

[©+©x@] y 2 

+ [©+©] F 3 
+[©] 01 + [@X®1 
— [®X@X@] 
+[®X©X®]=y"' (l) 

^ T n l 2n+i 
(2n + 4)! 

^ Tnil-k) 2 ^ 
^ (2n+4)! 

[~®x@]y 2 +[-®x®]F 3 
+[-©x®+©i 0i 

-[©X@X®1 

+[©x@x®i 

+[-®x®i=0 , (o 

1.6846-0.25672i 

(11.8000 — 3. 6695i) Y 2 
+ (10.2970-2.8566i) F 3 
+ (-5.4711 + O.93233i)0i 

10.807-0.12747i 

1.8198+0.004761 

(0.17030~0.66134i) Y 2 
+ (-0.09077-0.593411) Y 3 
+ (-0.41138+0.28864^)0, 

1.8685— 0.32229i 

(15.5930- 5.0935i)y 3 
+ (13.0080-3.7946i) F 3 
+ (-7.1158+ 1.3988i) 0i 

10.853-0.161401 

1.8110+0.00587t 

(-0.04177-0.87098*) Y 2 
+ (-0.23526 - 0.759481) Y 3 
+ (-O.514O3+O.37Ol7i)0i 
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FLUTTER OF A UNIFORM WING WITH AN ARBITRARILY PLACED MASS 
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Posjtive directions of axes and angles (forces and moments) are shown by arrows 


Axis A / . . 

/' /. 

Force; 
(parallel 
to axis) 
symbol 

J 

’"T i 1 

Moment about axis 

; ; Angle ’ 

)/ Velocities , 

! ' J*~ 

/y' 

, Designation 

Sym- 

bol 

Designation 

Syih- i 
bol 1 

Positive 

direction 

Designa- 
, tion 

“t • v A V 

Sym- 

bol 

Linear 
(cdmpo- 
nent along 
- ■ axis). 

• ■ *1 

Angular 

P '> 

Longitudinal-:.- 

■Lateral-.A ... 

-N ormaL— ... 

- *■ ' «‘ ; 

X 

Y 

: . Z 
\ ? 

X A 
Y 

. z 

Rolling— i— 
•' Pitching- A 
Yawing.-. 

' L 
v M 

' N : 

Y—^Z 
Z—^X 
’ X ->Y 

, ! 

Roll.-—. 

Pitch....... 

Yaw -A 

, A . • a \ v 

4> 

■9 

a:- 

■ • / 

' 

V 

w : 

r _<] 

V' • 

q 

T 


-Absolute coefficients of moment > ,, , 

, -• "Vy | L 1 Afy ■ ^7 ' iV 

, • . • ; Ll ~gps Cr ^‘acS , ~- Jn ~~qbS 

(rolling) ' (pitching) 1 (yawing) 


Angle of set of control surface (relative tb neutral 
position), 8. (Indicate surface by proper subscript.) 


"i'D-i Dianletcr " 
p Geometric pitch -y 

pJD Pitch' ratio A 
_ V r 'Inflow velocity ^ A . 

F,, Slipstream velocity 

T Thrust, absolute coefficient £7 


4. PROPELLER SYMBOLS 


P 


P Power, absolute coefficient C p7l ij Ji 


A P V 5 


: , t 


C, Speed-power coefficient — -’y 

7] A 1 Efficiency 

n Revolutions per second, rps 


v 




- pn 2 !}* t r % . t \ r ; ■ .A" • , A 1 / r ■- _/>y \ 

„ ■ o t ; / 4'- . - Effective helix angle = tarn H.wyzr I 

Q " ■; Torque, absolute coefficient A ■' ; . . \2ir rn/ j 

• ‘V _ s. numerical relations . ' 

1 lip— 76.04 kg -m/s— 550 ft-Ib/sec > . , 7 

I metric horsepower— 0.9863 hp 
1 mpli= 0.4470 mps 
.1 mps— 2.2369 mph ‘ ■ , 


\ . 


1 lb— 0.4536 kg; ' 

1 kg— 2.2046 lb ... - 

1 mi= 1,609.35 m-5,280 ft 
l.m= 3.2808 ft -• ■ 



